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ABSTRACT 

This article aims at proving the feasibility of the forecast of all the most relevant classi¬ 
cal atmospherical parameters for astronomical applications (wind speed and direction, 
temperature) above the ESO ground-base site of Cerro Paranal with a mesoscale at¬ 
mospherical model called Meso-Nh. In a precedent paper we have preliminarily treated 
the model performances obtained in reconstructing some key atmospherical parame¬ 
ters in the surface layer 0-30 m studying the bias and the RMSE on a statistical sample 
of 20 nights. Results were very encouraging and it appeared therefore mandatory to 
confirm such a good result on a much richer statistical sample. In this paper, the study 
was extended to a total sample of 129 nights between 2007 and 2011 distributed in 
different parts of the solar year. This large sample made our analysis more robust 
and definitive in terms of the model performances and permitted us to confirm the 
excellent performances of the model. Besides, we present an independent analysis of 
the model performances using the method of the contingency tables. Such a method 
permitted us to provide complementary key informations with respect to the bias and 
the RMSE particularly useful for an operational implementation of a forecast system. 

Key words: turbulence - site testing - atmospheric effects - methods: data analysis 
- methods: numerical 


1 INTRODUCTION 

This paper is part of a general study about the feasibility 
of the forecast of meteorological parameters and optical 
turbulence at ESO sites (Cerro Paranal and Cerro Arma- 
zones) in the framework of the MOSE project (Modeling 
ESO Sites). The MOSE project is pre s ented extensively 
in a previous paper dMasciadri et al.l l2013l l. We only 
recall here that the MOSE project aims at proving the 
feasibility of the forecast of all the most relevant classical 
atmospherical parameters for astronomical applications 
(wind speed intensity and direction, temperature, relative 
humidity) and the optical turbulence OT (C% profiles) 
with the integrated astro-climatic parameters derived from 
the Cff i.e. the seeing (e), the isoplanatic angle ( 9g ), the 
wavefront coherence time (to) above the two ESO sites of 
Cerro Paranal (site of the Very Large Telescope - VLT) and 
Cerro Armazones (site selected for the European Extremely 
Large Telescope - E-ELT). 

The final outcome of the project is to investigate the 
opportunity to implement an automatic system for the 
forecast of these parameters at the VLT Observatory at 
Cerro Paranal and at the E-ELT Observatory at Cerro 
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Armazones._ 

In a previous paper dLascaux et al. IB), we presented 
results of a deep analysis of the bias and RMSE between 
observations and model outputs of absolute temperature, 
wind speed and direction above both astronomical sites 
(Paranal and Armazones). These statistical operators 
provide fundamental informations on systematic and sta¬ 
tistical model errors. This statistical study was performed 
on a sample of 20 nights, the same for each site. A part 
of the nights is co ming from the PAR2007 campaign 
(IDali Ali et al. 120101 1. and the sample was completed with 
nights with available observations from the same period of 
the year. In the same study we also provided a detailed 
analysis on the model performances on individual nights. 
The most important results we obtained can be summarized 
as in the following. The model showed a very good score of 
success in terms of bias and RMSE for absolute temperature 
(median bias in the [0.03,0.64]°C range and median RMSE 
in the [0.64,0.93]°C range) and wind direction (median bias 
in the [-1.01,-8.55]° range and median RMSE in the [30,41]° 
range for a wind speed threshold of 2 m s^ 1 equivalent to 
a RMSEj-eiatiue in the [17,23]% range) when it is used in 
grid-nesting configuration and the horizontal resolution of 
the innermost model domain is 500 m. To obtain equivalent 
satisfactory results for the wind speed we need to use a 
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horizontal resolution of the innermost model domain equal 
to 100 m. Under this condition we obtained a median bias 
within 0.93 m s -1 and a median RMSE within 2.18 m s^ 1 . 
These results were, in conclusion, very promising. 

In th is p aper the preliminary analysis done by 
lLascaux et al.l ( 2013 1 was extended to a sample of 129 
nights uniformly distributed between 2007 and 2011. The 
analysis is performed for the Cerro Paranal case. It was 
possible to have access to homogeneous measurements 
distributed over such a long period of time, only at this 
sit^B All the 129 nights belong to different periods of 
the year, and not onl y to summ er (like for the sample 
of 20 nights from lLascaux et alj (l2013ll 1. This permits 
us to analyze definitely the robustness of the model on 
a rich statistical sample, independently from the period 
of the year considered. All the measurements, at Cerro 
Paranal, are part of t he VLT Astronomical Site Monitor 
dSandrock fc Amesticall2009l) . 


is definitely very low (typically less than 5 nights for year). 
In all the other nights the RH does not really represent a 
critical issue in terms of scheduling of observations because 
the RH value is not high enough to affect the observations. 
The study of the relative humidity close to the ground is 
therefore not very critical and important for telescopes in 
this region. Also it has been observed in occasion of the 
TMT site characterization in this region (T. Travouillon, 
private communication), that the reliability of the relative 
humidity measurements, when values are so low, can hardly 
be assured. This makes difficult (and mostly useless) to 
investigate the model performances in reconstructing the 
RH on the same sample of nights selected for the analysis 
of the other parameters. A more interesting analysis might 
be to check the model performances in reconstructing the 
RH values in those few nights in each year where RH is 
higher than the threshold imposing to close the dome of the 
telescope. In Annex we report satisfactory results obtained 
in this sense. 


Besides that, on the same sample of nights, we present 
a complementary statistical analysis based on a different 
approach. Indeed, it is worth to highlight that bias and 
RMSE, in spite to be fundamental statistical operators 
providing key informations on model performances, do not 
answer to the necessary information one wou ld lik e to have 
in terms of model performances |T hornes & Stenhensonl 
l200ll ; lNurmill2003l : I Jolliffe fc Stephensonl 2003f l. A key rule 
in the complex art of the model prediction estimation is 
that it makes no sense to quantify the model performance 
with one parameter only or one method only. To investigate 
the quality of a model prediction, a method widely used in 
physics of the atmosphere as well as in other fields such 
as economy and medicine, consists in constructing and 
analyzing contingency tables. From these tables one can 
derive a number of different parameters that describe the 
quality of the model performances. 

A contingency table allows for the analysis of the rela¬ 
tionship between two or more categorical variables. One 
of the first classic examples of a forecast verification usin g 
contingency tables is the tornado forecast from (Finley 
). A contingency table is a table with nxn entries 


that displays the distribution of modeled outputs and 
observations in terms of frequencies or relative frequencies. 
Here the variables we considered are the observations 
of different categories of temperature, wind speed, wind 
direction, near the surface, and the respective reconstructed 
parameters by the Meso-NH model. 


In this part of Chile the weather is particularly dry, 
the relative humidity (RH) in the surface layer is typically 
below 30%. The frequency of the nights with a relative 
humidity above 80%, level at which the dome of the 
telescope has to be closed and observations are not allowed, 


1 At Cerro Armazones during this period part of measurements 
have been done with TMT instrumentation and part with ESO 
one with sensors located at a different number of levels and 
heights above the ground. This lack of homogeneity in measure¬ 
ments did not permit us to perform the analysis at Armazones 
using the same procedure. We therefore have decided in this paper 
to focus our attention only on Cerro Paranal. 


In Section0we describe the model configuration defined 
for this study. In Section [3] we report the bias and RMSE 
obtained for the new sample of 129 nights, at Cerro Paranal. 
In Section [I] we define what is a contingency table for dif¬ 
ferent discretized values n of observations and predictions 
and which parameters qualifying the model behavior can be 
retrieved from them. In Section [S] we describe the criteria 
we used to identify the thresholds separating the n cate¬ 
gories. In Sections [b] [ 7 ] and [8] we provide the results in terms 
of model performances (reconstruction of the absolute tem¬ 
perature, the wind speed and direction in the 0-30 m range). 
Finally conclusions are drawn in Section [9] 


2 BRIEF OVERVIEW OF THE NUMERICAL 
CONFIGURATION 

2.1 Sample selection 

Th e 129 ni ghts ha ve bee n selected starting with the 20 nights 
of lLascaux et al] J2013ll . Those 20 nights were all concen¬ 
trated i n Novem ber and December 2007 (the reader can refer 
to lLascaux et al] (l2013l l to understand how this first sample 
was constructed). When results obtained with the sample 
of 20 nights indicated a very promising model behavior, to 
have a better evaluation of the model performances in sta¬ 
tistical terms, it was decided to extend this sample to a 
more consistent one. Two different periods, distant in time, 
were arbitrarily chosen (the only constraint was the avail¬ 
ability of continuous measurements): 1/from January 2007 
to December 2007 and 2/from June 2010 to May 2011. We 
tried to select nights that were regularly spaced in time, 
as much as possible, in order to cover all the four seasons 
(which are characterized by slightly different meteorological 
behaviours). The result was 55 nights in 2007 (the 20 ones 
mentioned above, plus around 3 nights per month), and 74 
nights in 2010-2011 (around 6 nights per month). No spe¬ 
cific criterium has been used on the goodness of the weather 
conditions. In other words, nights include good and/or bad 
weather conditions therefore no specific bias affects the sam¬ 
ple. All these nights were characterized by different condi¬ 
tions of wind speed, temperature and relative humidity (see 
discussion in Section ®. 
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Table 1. Meso-NH model configurations. In the second column the horizontal resolution AX , in the third column the number of grid 
points and in the fourth column the horizontal surface covered by the model domain. 


Domain 

AX (km) 

Grid Points 

Domain size 
(km) 

Domain 1 

10 

80x80 

800x800 

Domain 2 

2.5 

64x64 

160x160 

Domain 3 

0.5 

150x100 

75x50 

Domain 4 

0.1 

100x100 

10x10 


2.2 Model configuration 

All the numerical simulations of the nights (129 in to¬ 
tal) presented in this study have been performed with the 
meso scale numerical weather model Meso- Nffl (iLafore et ali 
11998h and the Astro-Meso- NH code for the optical turbu¬ 
lence (iMasciadri et all 19991 '). The model has been developed 
by the Centre National des Recherches Meteorologiques 
(CNRM) and Laboratoire d’Aerologie (LA) from Universite 
Paul Sabatier (Toulouse). The Meso-NH model can simulate 
the temporal evolution of three-dimensional meteorological 
parameters over a selected finite ar ea of the globe. We refer 
the reader to lMasciadri et all (l2013lh Sec. 3.3, for the general 
model configuration and the physical packages used for this 
study. We just recall her e that we used the grid-nesting tech¬ 
nique dStein et al.ll20odh . that consists in using different im¬ 
bricated domains of the Digital Elevation Models (DEM i.e 
orography) extended on smaller and smaller surfaces, with 
increasing horizontal resolution but with the same vertical 
grid. In this study we use two different configurations. The 
first grid-nesting configuration employed three domains (Ta¬ 
ble [l]) and the innermost resolution is AA = 500 m. The 
second configuration is made of four imbricated domains, 
the first same three as the previous configuration, plus one 
centered at Cerro Paranal site, with a horizontal resolution 
of AA' = 100 m (all domains of Table |TJ . Along the z-axis 
we have 62 levels distributed as follows: a first vertical grid 
point equal to 5 m, a logarithmic stretching of 20 % up to 
3.5 km above the ground, and an almost constant vertical 
grid size of ~600 m up to 23.8 km. 

All the simulations done for the analyses discussed in this 
paper, were initialized the day before at 18 UT and forced 
every 6 hours with the analyses from the ECMWF (Eu¬ 
ropean Centre for Medium-Range Weather Forecasts), and 
finished at 09 UT (05 LT) of the simulated day (for a total 
durati on of 15 hours). The reader can refer to lLascaux et al.l 
(l2013h to have more information about the site coordinates 
and ground altitude. 

Before computing bias and RMSE, and before construct¬ 
ing the contingency tables, all the investigated meteorologi¬ 
cal parameters (measurements and Meso-NH outputs) have 
been first averaged with a moving average of 1 hour to cut 
the high frequencies and then resampled over a 20-minute 
interval. The moving average of one hour has been selected 
because astronomers are in reality more interested in the 
trend of the prediction. It has been observed that the moving 
average is more efficient to identify the model and measure¬ 
ments trends than a simple re-sampling. The high frequency 


variability on a time scale of 5 or 10 minutes is less relevant 
and useless to predict because the scheduling of scientific 
programs could not be tuned with such high frequency. As¬ 
tronomers are interested in identifying if the trend of a pa¬ 
rameter is increasing, decreasing or is stationary to be able 
to take a decision in changing a modality of observation or 
a scientific program. The selection of the interval of 20 min¬ 
utes for the re-sampling of measurements is justified by the 
fact that this is more or less the effective time necessary to 
switch from a modality of observation to another one. The 
values of one hour and 20 minutes have been selected in 
agreement with ESO staff. In the analysis presented in this 
paper we will consider the model abilities in reconstructing 
the atmospherical parameters during nighttime onljjf). 


3 STATISTICAL ANALYSIS 


In this section we estimate the statistical model reliability 
using three statistical operators: the bias, the root mean 
square error (RMSE) and the bias-corrected RMSE. 


N 

BIAS = J2 

i= 1 


(AQ 

N 


(1) 


RMSE = 



( A ,) 2 

N 


( 2 ) 


with Aj = Yi—Xi where X\ are the individual observations, 
Yi the individual simulations parameters calculated at the 
same time and N is the number of times for which a couple 
(A i,Yi) is available with both A \ and Yi different from zero. 
Because the wind direction is a circular variable, we define 
A i for the wind direction as: 


j Yi - Xi 

A, = I Yi - AL - 360° 


if | Yt -Xi\< 180° 
if Yi - Xi > 180° 


Yi — Xi + 360° if Yj - Xi < -180° 


(3) 


From the bias and the RMSE we deduce the bias-corrected 
RMSE (it): 

a = sjRMSE 2 - BIAS 2 (4) 


which represents the error of the model once the systematic 
error (the bias) is removed. For the wind direction we add 


2 http://mesonh.aero.obs-mip.fr/mesonh/ 


3 To perform daytime predictions, a different model configuration 
should be selected. 
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Figure 1. Scattered plots of the absolute temperature (top row), of the wind speed with the AX = 500 m configuration (second row), 
of the wind speed with the AX = 100 m configuration (third row) and of the wind direction (bottom row). A moving average (1 hour) 
have been applied on both measurements and Meso-NH outputs (cf. text). Each point represents an average of 20 minutes. The thin 
black line is the regression line passing by the origin. The red dots correspond to points for which the observed wind was inferior to 
WSth = 3 m-s —1 . The corresponding values of bias and RMSE of the wind direction for the sample without these dots, are reported in 
red in the bottom right of the wind direction scattered plots. The values of the biases and the RMSE are reported in Table [2] 
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Table 2. Near surface bias, RMSE, bias-corrected RMSE a (Meson-NH minus Observations), at Cerro Paranal, for the 129 nights sample 
(cf. Figure [l]). The relative RMSE is also reported in the case of the wind direction. The computations have been made considering the 
whole sample of data points. For the wind direction, data corresponding to a wind speed inferior to 3 m-s —1 have been discarded from 
the computations. 


CERRO 

PARANAL 

Temperature (°C) 
AX = 500 m 

Wind speed (m-s 1 ) 
AX = 500 m 

Wind speed (m-s 1 ) 
AX = 100 m 

Wind direction (°) 
AX = 500 m 

2 m 

30 m 

10 m 

30 m 

10 m 

30 m 

10 m 

30 m 

BIAS 

-0.12 

-0.37 

-2.90 

-1.39 

-1.16 

-0.24 

3.4 

-0.2 

RMSE 

0.99 

1.02 

4.03 

3.04 

2.73 

2.70 

36.7 

35.7 

17 

0.98 

0.95 

2.80 

2.70 

2.47 

2.69 

36.5 

35.7 

RMSE re / 







20.4% 

19.8% 


Table 3. Individual nights model performance. Near surface median bias, RMSE and bias-corrected RMSE (Meson-NH minus Obser¬ 
vations). In small fonts, the 1st and 3rd quartiles. The values for the wind speed were obtained with the AX = 100 m configuration for 
Meso-NH, whereas the values obtained for the temperature and the wind direction were obtained with AX = 500 m configuration. 
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Table 4. Winter: near surface median bias, RMSE and bias-corrected RMSE (Meson-NH minus Observations). Only the 54 winter nights 
are considered. In small fonts, the 1st and 3rd quartiles. The values for the wind speed were obtained with the AX = 100 m configuration 
for Meso-NH, whereas the values obtained for the temperature and the wind direction were obtained with AX = 500 m configuration. 
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Table 5. Summer: near surface median bias, RMSE and bias-corrected RMSE (Meson-NH minus Observations). Only the 75 summer 
nights are considered. In small fonts, the 1st and 3rd quartiles. The values for the wind speed were obtained with the AX = 100 m 
configuration for Meso-NH, whereas the values obtained for the temperature and the wind direction were obtained with AX = 500 m 
configuration. 
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another operator, the relative RMSE: 

RMSErel = x 100% (5) 

lot) 

with RMSE expressed in degrees. This is justified because 
the wind direction is a circular variable and the maximum 
error is obtained when RMSE is equal to 180°. 

Figure [T] shows the scattered plots for the temperature, the 
wind speed (with both AX = 500 m and AX = 100 m config¬ 
urations) and the wind direction (with and without filtering 
the lowest wind speed). The first thing to notice is the total 
extent of the observed values for all the meteorological pa¬ 
rameters. The temperature goes from around 0°C to 20°C, 
the wind speed reaches 20 m-s -1 and the wind direction 
covers the whole [0-360]° range. This is in agreement with 
the minimum and the maximum abserved values during a 
whole year at the site, as is reported on the Astoclimatol- 
ogy Website of ParanaJj. This assures us that all the possible 
different meteorological conditions were encountered for the 
129 simulated nights. 

The forecast of the absolute temperature and the wind 
direction near the surface at Cerro Paranal, with the Meso- 
NH model, are excellent. Considering the whole sample (Fig- 
ure[7]and Table[2|, the bias in the temperature is very small 
(not larger than 0.37°C in absolute value) and the RMSE is 
not larger than 1°C. 

The bias in the wind direction is almost zero, and the 
RMSE is around 37° (the data for which the wind velocity 
was inferior to 3 m-s -1 have been excluded from the com¬ 
putations. These data are associated in general to a great 
dispersion for the wind direction but are not of interesting 
data for astronomers). This corresponds to a RMSE re ; of 
around 20%. 

For t he wind speed , we o bserve that the conclusions 
found in lLascaux et alJ ( 2013 : are confirmed. The high 
horizontal resolution (AX = 100 m) is mandatory to reduce 
the bias to acceptable values (from 2.90 m-s -1 to 1.16 m-s -1 
at 10 m and from 1.39 m-s -1 to 0.24 m-s -1 at 30 m, in 
absolute values). With the AX = 100 m resolution, the 
RMSE is around 2.7 m-s -1 . 

When looking at the single nights statistics, the per¬ 
formances are even better. We computed the bias, RMSE 
and bias-corrected RMSE for every night of the sample, 
and constructed the cumulative distributions at every levels 
for every meteorological parameter. From the cumulative 
distributions we can extract the median and the first and 
third quartiles. Figure [ATI and lA2l in the appendix show the 
cumulative distributions at all levels. Tabic 0 summarizes 
all the median values, together with the first and third 
quartiles, for the single nights bias, RMSE and a, for the 
absolute temperature, the wind speed with the AX = 100 m 
configuration, and the wind direction. The median RMSE 
for the absolute temperature is less than 0.92°C. It is 
inferior to 2.30 m-s -1 for the wind speed, and inferior to 
29° for the wind direction. 

To study the seasonal variation and check if some 
systematic better/worst model performance is observed 
in summer and winter we also calculated bias, RMSE 


and bias-corrected RMSE for the two periods of [April- 
September] that we call winter and [October-March] that 
we call summer (Table []] and Table [5j). The two periods 
are made of 54 and 75 nights, respectively. Summer period 
has more nights bec ause of the p re sence of the first 20 
nights analyzed by lLascaux et al.l (l2013f) that were all 
concentrated in the months of November and December 
2007. We observe that, for the temperature, in summer 
we have a slightly better median bias but a slightly larger 
median RMSE as well as sigma (bias-corrected RMSE) than 
in winter. In any case, the median RMSE remains below 
1°C. For the wind direction the median bias of the absolute 
temperature is sightly better in summer than in winter 
but the median RMSE is larger (almost the double 33-35° 
instead of 16-19°) in summer than in winter. The sigma 
also shows the same trend even if the difference between 
summer and winter is less marked. For the wind speed negli¬ 
gible difference are observed for the median bias and RMSE. 

Now that we have characterized the performances of the 
model using standards statistical parameters (bias, RMSE 
and a), we want to know how good or bad is the model 
in predicting some specific events (like strong wind speed) 
or categories (like the quadrants from which the wind is 
blowing). For this we have constructed contingency tables. 


4 http://www.eso.org/gen-fac/pubs/astclim/paranal 
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4 WHAT IS A CONTINGENCY TABLE? 

Like mentioned in the Introduction, a contingency table al¬ 
lows for the analysis of the relationship between two or more 
categorical variables. Table [6] is an example of a generic 2x2 
contingency table, where the observations are divided in 2 
categories (YES and NO), that could for example corre¬ 
spond to a wind above a given threshold (YES) and a wind 
below this threshold (NO). In the 2x2 contingency table, 
four combinations are possible: 

-’a’ is the number of times an event was well reproduced by 
the model (hit); 

-’6’ is the number of times an event was erroneously fore¬ 
casted by the model (false alarms); 

-V is the number of times an event was not reproduced by 
the model (miss); 

-’d’ is the number of times the absence of an event was well 
reproduced by the model (correct rejection). 

Using a, b , c and d (and N = a + b +c + d), we can compute 
different probabilities useful to have an insight on how well 
(or bad) the model performed for a particular event. 

Hereafter we list all the different simple scores we will 
use in the following of the article (to define them we use the 
numbers a, b, c, d and N from the generic 2x2 contingency 
table of Table [5]): 

-the percent of correct detections PC (in %): 

PC = x 100; 0% < PC < 100% (6) 

PC =100% is the best score, and corresponds to a perfect 
forecast. The PC alone is not sufficient because sometimes 
it can be influenced by the most common category and it 
does not provide us informations on the ability of the model 
in detecting observations in the different categories. It is 
therefore always preferable to provide also the probability 
of detection ( POD ) together to the PC. 

-the probability of detection (POD, in %) of a given event, 
or hit-rate: 

POD(event i) = —— x 100; 0% < POD < 100% (7) 

a + c 

POD =100% is the best score. POD measures the propor¬ 
tion of observed events that have been correctly predicted 
by the model. In the 2x2 contingency table, only one event 
is generally considered, and thus only one POD is calcu¬ 
lated. However we could also consider the non occurrence of 
the event as an event itself, and then define a second POD 
as: 

POD(eventi ) = x 100; 0% < POD < 100% (8) 

It can be of interest to highlight that, for a total ran¬ 
dom prediction and in the case of a 2x2 contingency table, 
a = b = c = d. That means that all POD are equal to 50%, 
and PC'=50%. The model is useful if it performs better than 
this random case. 

A huge number of others statistical parameters can be de¬ 
duced from a contingency table (such as false alarm ratio, 
probability of false detection, etc.). We will limit the current 
study to the analysis of the parameters mentioned above 
that already provide a consistent panorama of the model 
potentialities. 

For our purpose a 2x2 table is, however, a too simplified 
analysis. A 3x3 table is definitely more appropriated. It 


consists in dividing the observed values in three categories 
delimited by some thresholds. An example of a 3x3 contin¬ 
gency tables is shown in Table [7] Equivalent PC and POD 
can be defined using a, b, c, d, e, f, g, h, i and N from 
Table [3 In the case of the 3x3 contingency tables, we also 
add a third parameter we call EBD (for extremely bad de¬ 
tection). 


pC = a + e + i xW0 

N 

(9) 

POD(eventi) =- - - x 100 

v ' a + d + g 

(10) 

POD(event 2 ) = ---—- x 100 

b + e + h 

(11) 

POD(event 3 ) =- - - x 100 

c + f +1 

(12) 

EBD = C + 5 x 100 

(13) 


EBD represents the percent of the most distant predictions, 
by the model, from the observations. For a perfect forecast, 
it is equal to 0%. In the case of a perfectly random forecast 
(a = b = ... = i = -j), all POD are equal to 33%, PC= 33%, 
and EBD= 22.2%. These values are a good reference to eval¬ 
uate the performances of the model. 

For the wind direction a 3x3 table is not appropriated be¬ 
cause this is a circular variable. A 4x4 (divided in four quad¬ 
rants with 90° each) is therefore a more suitable solution (see 
Table [8]). PC and POD are defined in this case as: 


PC = 100x a + f + k+P 

N 

(14) 

POD(event \) = 100 x --- 

a + e + i + m 

(15) 

f 

POD(event 2 ) = 100 x ---- : - 

V ’ b+f+j+n 

(16) 

POD(event7) = 100 x - — - 

c+g+k+o 

(17) 

POD(eventi) = 100 x —— P - 

d + h + l + p 

(18) 

EBD = 100x C + h + i + n 

N 

(19) 

If we consider the random case (a = b = ... = p 

= K) 
16'’ 


jljl vv*3 LUiiniuci Lite jl vaoc yu, — u — ... — jj — yg 

PC= 25%, all POD are equal to 25%, and EBD= 25%. 


From this point and until the end of the article, we will write 
PODi instead of POD(eventi) with i the event considered. 


5 CLIMATOLOGICAL TERTILES AT CERRO 
PARANAL 

To construct the contingency tables that we will discuss in 
the next sections, it is necessary to define some thresholds 
limiting the discretized intervals of the observed values and 
to divide the sample in three categories. In the cases of 
the 3x3 contingency tables, we need two thresholds. In this 
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Table 6. Generic 2x2 contingency table. 

EVENT 


jJ YES 

H 

Q 

O 

2 NO 


OBSERVATIONS 


YES 

NO 

Total 

a 

b 

a+b 

(hit) 

(false alarm) 

Yes (Model) 

c 

d 

c+d 

(miss) 

(correct rejection) 

No (Model) 


Total 


a+c 

Yes (OBS) 


b+d 

No (OBS) 


N=a+b+c+d 
Total of events 


Table 7. Generic 3x3 contingency table. 

Intervals 


OBSERVATIONS 

1 2 3 


Total 


H 

D 

O 


a 

(hit 1) 


(hit 2) 
h 


(hit 3) 


a+b+c 

1 (Model) 

d+e+f 

2 (Model) 


g+h+i 
3 (Model) 


, a+d+g b+e+h c+f+i N=a+b+c+d+e+f+g+h+i 
° a 1 (OBS) 2 (OBS) 3 (OBS) Total of events 


Table 9. Climatological tertiles for the wind speed at Cerro 
Paranal (left column: 33%, central column: median, right column: 
66 %) 


Cerro 

Wind Speed (in m s 

- 1 ) 

Paranal 

33% 

Median (50%) 

66% 

10 m 

3.98 

5.47 

7.25 

30 m 

4.40 

6.10 

8.10 


Table 10. Climatological tertiles for the absolute temperature at 
Cerro Paranal (left column: 33%, central column: median, right 
column: 66%) 


Cerro 

Temperature (in Celsius) 

Paranal 

33% 

Median (50%) 

66% 

2 m 

11.03 

12.24 

13.27 

30 m 

11.40 

12.58 

13.64 


study, we decided to use the climatological tertiles computed 
thanks to the available measurements. We used 6 complete 
years of measurements from a mast, from January 1, 2006 to 
December 31, 2012, that are part of the VLT Astronomical 
Site Monitor (ISandrock fe Amestical 120091) . At each level, 
and for each parameter, we computed the median value, 
and the tertiles, of all the available measurements. Table [9] 
shows the results for the wind speed at Cerro Paranal. Ta¬ 
bles ED shows the results for the absolute temperature at 
Cerro Paranal. 


6 ABSOLUTE TEMPERATURE 

The first meteorological parameter reconstructed by the 
model that we investigate is the absolute temperature near 
the ground. We report the results only for the standard 
model (AX=500 m). We highlight that results obtained 
with the high horizontal resolution model (AX=100 m) (not 


shown here) are very similar. Tabic fill and Table fl2l are the 
contingency tables for the absolute temperature at Cerro 
Paranal, at 2 m and 30 m, respectively. Each estimate cor¬ 
responds to the average over a 20-minutes interval (as indi¬ 
cated in Section E3>. For sake of simplification, we used for 
thresholds, the rounded values taken from the climatoligical 
study of Section [5] We can observe that the percent of good 
forecast by the model (characterized by PC) is very good 
at every level. PC is between 73% and 75%. More over, an 
EBD very close to 0% at both levels, is the sign that the 
model never produces extremely bad forecasts for the tem¬ 
perature. If we look at the probability of detection of the 
temperature in a given interval (characterized by POD), 
once again the model demonstrates its accuracy. POD is 
between 55.2% (detection of temperature between 11.5°C 
and 13.5°C at 30 m) and 93.7% (detection of a temperature 
inferior to 11.5°C at 30 m). In all cases, PC and PODs are 
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Table 8 . Generic 4x4 contingency table. 


Intervals 


OBSERVATIONS 

12 3 


4 


Total 


a 

(hit 1) 


d 


a+b+c+d 
1 (Model) 


O 

s 


2 


3 


e 


f 

(hit 2) 


h 


k 

(hit 3) 


1 


e+f+g+h 

2 (Model) 

i+j+k+1 

3 (Model) 


p m+n+o+p 

(hit 4) 4 (Model) 


Total a+e+i+m b+f+j+n c+g+k+o d+h+l+p 
1 (OBS) 2 (OBS) 3 (OBS) 4 (OBS) 


N=a+b+c+d+e+f+g+h 
+i+j+k+l+m+n+o+p 
Total of events 


Table 11. 3 x3 contingency table for the absolute temperature during the night, at 2 m a.g.l. at Cerro Paranal, for the sample of 129 
nights. We use the Meso-NH AX = 500 m configuration. 


Division by tertiles (climatology) OBSERVATIONS 

C. Paranal-2 m T < 11°C 11°C <T < 13.5°C T > 13.5°C 


T < 11°C 

1157 

315 

0 

11°C < T < 13.5°C 

96 

972 

258 

T > 13.5°C 

8 

178 

499 


Total points = 3483; PC= 75.5%; EBD= 0.2% 
PODi=91.8%; POD 2 = 66.3%; POD 3 = 65.9% 


well larger than 33% (value of random case). This proves 
therefore the utility of the model. 


7 WIND SPEED 


Considering the tendency of the model to slightly underes¬ 
timate the wind velocity (see section [3]), (even if this un¬ 
derestimation is obviously quantitatively different between 
AX=500 m and AX=100 m), it was decided to evaluate the 
predictive power of the model using a forecasted wind speed 
corrected by the multiplicative bias Bmuit computed as: 


Bmult (^) 


WS mod (k) 

WsFF(k) 


( 20 ) 


with WSm 0 d(k) the average of the forecasted wind speed 
over the whole sample of the 129 nights, and WS 0 b a (fc) the 
average of the observed wind speed from the same sample, 
at the level k (10 m and 30 m for Cerro Paranal). The dif¬ 
ferent values are summarized in Table I~HT1 The choice of the 
multiplicative bias is justified by the fact that we want the 
highest values of the wind speed to be more corrected than 
the lowest values. It is important to note that the values of 
Bmuit(fc) can be considered as constant, whatever the sub¬ 
sample considered. Indeed, we have computed it with the 


Table 13. Wind speed multiplicative biases (from Eg. 1201 . Values 
are rounded. 




10 m 

30 m 

C. 

Paranal 

0.5 

0.8 

AA' 

= 500 m 

C. 

Paranal 

0.8 

0.9 

AA' 

= 100 m 


initial sample of 20 nights, with another subsample of 73 
nights, and with the final sample of 129 nights, and in all 
cases its values were the one reported in Table m (varia¬ 
tions of the second digit only are noticeable, and thus not 
reported here). This tells us that is useless to consider a 
calibration sub-sample and a complementary and indepen¬ 
dent validation or testing sub-sample. It is therefore more 
appropriate to apply this multiplicative factor on the whole 
sample night. 
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Table 12. 3x3 contingency table for the absolute temperature during the night, at 30 m a.g.l. at Cerro Paranal, for the sample of 129 
nights. We use the Meso-NH AX = 500 m configuration. 

Division by tertiles (climatology) OBSERVATIONS 

C. Paranal - 30 m T < 11.5°C 11.5°C < T < 13.5°C T > 13.5°C 


T < 11.5°C 

1250 

422 

2 

11.5°C < T < 13.5°C 

68 

663 

292 

T > 13.5°C 

16 

117 

653 


Total points = 3483; PC=73.7%; EBD= 0.5% 
PODi=93.7%; POD 2 = 55.2%; POD 3 =6 9.0% 


POD vs Vlim 



Figure 2. Evolution of the probability of detection (POD) from 
2x2 contingency tables, by Meso-NH, of winds stronger than 
Vlim > i n function of Vum- Dashed lines, with AX = 100 m con¬ 
figuration. Continuous lines, with AX = 500 m configuration. In 
all cases, the Meso-NH wind speed was corrected by the multi¬ 
plicative bias of Table na CP stands for Cerro Paranal. 

7.1 Sensibility of the model detection to the wind 
speed threshold 

Considering that the most critical issue for the wind 
speed is the model performances in detecting the strong 
wind speed we performed the following exercise before to 
calculate the 3x3 contingency tables. To see how the model 
outputs depend on the choice of the threshold used for the 
wind speed, we plotted on Fig. [2] the dependence of POD\ 
from 2x2 contingency tables, on the threshold Vum, for 
both AX = 500 m and AX = 100 m configurations. Vum 
varies between 4 m-s -1 and 11 m-s -1 . The two Meso-NH 
configurations (AX = 500 and 100 meters) were used with 
the correction of the wind by the multiplicative bias. The 
first thing we can notice in Fig [2] is that, the performance 
of the simulations with the AX = 100 m configuration 
depends less on the threshold than the performance of 
the AX = 500 m configuration. Up to Vu m = 10 m-s _1 , 
with AX = 100 m all POD\ are very good with the 
worst PODs at 10 m-s -1 between 70% and 80%. On the 
contrary, the model performance decreases faster with Vu m , 
for the AX' = 500 m configuration. At Vu m = 9 m-s _1 , 
the performance is still satisfactory (all POD i are around 
60%), but for Vum > 9 m-s -1 , the performances are 


Table 14. Contingency table considering the event: winds 
stronger than 6 m-s— U The site is Cerro Paranal. The altitude 
is 10 m a. g. 1. The wind speed is computed with the Meso-NH 
AX = 100 m configuration and corrected by the multiplicative 
bias from Table fT3l 

WS > 6 m-s" 1 OBSERVATIONS 

C. Paranal - 10m YES NO 


A 




a 

YES 

1230 

418 

Q 




u 

§ 

NO 

363 

1445 


Total points = 3456 
PC=77A% 

PODi(WS>6m-s- 1 )=77.2% 

POD 2 (WS<6m-s- 1 )=77.6% 

degraded (all POD\ are inferior to 50%, down to ~30% 
for Vum = 11 m-s -1 ). This tells us that the 100 m con¬ 
figuration guarantees a better model performances than 
the 500 m one. As an example, if we consider at Cerro 
Paranal a reasonable value of Va m =6 m-s -1 and we study 
the probability to detect wind speed stronger than this 
threshold we obtain very satisfactory results (Table fill) for 
the PODs that are typically of the order of 77%. 

The interest for astronomers is mainly to be able to 
identify conditions of strong wind and this result tells us 
that such a model in this configuration well answers to this 
necessity. 

7.2 Wind speed 3x3 contingency tables 

In this section we investigate the ability of the model in fore¬ 
casting the wind speed as we did in Section [6] As we have 
seen in Section rm we consider the best configuration for 
the model i.e. the high horizontal (AX=100 m) resolution 
configuration. The contingency tables have been divided in 
three intervals, for each level concerned, using the rounded 
climatological tertiles from section [5] as thresholds. They are 
reported in Tables [Tsl and flGl at 10 m and 30 m, respectively. 
PC is significatively better than 33% (random case): it is 
60.0% at 10 m and 60.6% at 30 m. The POD also are very 
satisfactory, between 47.5% and 79.0%. More over, the EBD 
are negligible (less than 4.4%). If we consider POD 3 only 
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Table 15. 3x3 contingency table for the wind speed during the night, at 10 m a.g.l. at Cerro Paranal. We use the Meso-NH AX = 100 m 
configuration with the wind corrected by the multiplicative bias. 

Division by tertiles (climatology) OBSERVATIONS 

C. Paranal - 10 m T < 4 m-s~ 1 4 m-s~ 1 < T < 7 m-s~ 1 T > 7 m-s ~ 1 


T < 4 m-s ~ 1 

641 

283 

50 

> 1 < T < 7 Tfl'S 1 

395 

510 

269 

T > 7 m-s ~ 1 

103 

281 

924 


Total points = 3456; PC=60.0%; EBD= 4.4% 
POP>i=56.2%; POD 2 = 47.5%; POD 3 = 74.3% 


Table 16. 3x3 contingency table for the wind speed during the night, at 30 m a.g.l. at Cerro Paranal. We use the Meso-NH AX = 100 m 
configuration with the wind corrected by the multiplicative bias. 

Division by tertiles (climatology) OBSERVATIONS 

C. Paranal - 30 m T < 4 m-s~ 1 4 m-s~ 1 < T < 8 m-s~ 1 T > 8 m-s ~ 1 


T < 4 m-s ~ 1 

554 

289 

27 

i 1 < T < 8 rri'S 1 

400 

559 

234 

T > 8 m-s ~ 1 

101 

311 

981 


Total points = 3456; PC=60.6%; EBD= 3.7% 
PODi=52.5%; POD 2 = 48.2%; POD 3 = 79.0% 

(that is, detection of the wind superior to a fixed thresh¬ 
old, the 2 nd climatological tertile), that is the most interest¬ 
ing informations for astronomers, the results are even better 
(-79.0%). 

A strong wind is indeed the main cause of vibrations of the 
primary mirror and adaptive secondary (a critical element 
of the adaptive optics system). All POD 3 for both sites and 
at all altitudes are excellent, between 74.3% and 79.0%. This 
demonstrates the ability of the model in predicting at Cerro 
Paranal critical wind speed values for AO applications. 


8 WIND DIRECTION 

In this section we investigate the reconstruction of the wind 
direction by the Meso-NH model with the AX = 500 m 
configuration. We have constructed 4x4 contingency tables 
divided in four quadrants of 90° each. 


N 

0 ° 



Figure 3. Convention chosen for the wind direction (0° < a < 
360°). The black quadrants correspond to the intervals of 90° cho¬ 
sen for the 4x4 contingency tables of Table [171 the red quadrants 
correspond to the intervals of 90° chosen for the 4x4 contingency 
tables of Table [l8l The hatched section is discussed in the text. 
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Table 17. 4 x4 contingency table for the wind direction a during 
the night, at 10 m, and 30 m a.g.l. at Cerro Paranal. We use the 
Meso-NH AX = 500 m configuration. We filter out the observed 
wind inferior to 3 NE corresponds to 0° < a < 90°; 

SE corresponds to 90° < a < 180°; SW corresponds to 180° < 
a < 270°; NW corresponds to 270° < a < 360°. a is the wind 
direction (see Fig. 


Table 18. 4x4 contingency table for the wind direction a during 
the night, at 10 m, and 30 m a.g.l. at Cerro Paranal. We use the 
Meso-NH AX = 500 m configuration. We filter out the observed 
wind inferior to 3 m-s —1 . N corresponds to —45° < a < 45°; E 
corresponds to 45° < a < 135°; S corresponds to 135° < a < 
225°; W corresponds to 225° < a < 315°. a is the wind direction 
(see Fig. 0. 


OBSERVATIONS OBSERVATIONS 


C. Paranal - all levels 

N-E 

S-E 

S-W 

N-W 

C. Paranal - all levels 

N 

E 

S 

W 

N-E 

2809 

212 

35 

815 

N 

3807 

536 

68 

138 

A 





A 





H S-E 

102 

734 

109 

28 

H E 

196 

473 

119 

4 

P 





Q 





g S-W 

23 

87 

63 

11 

S 

18 

141 

512 

25 

N-W 

505 

87 

19 

614 

W 

83 

30 

38 

64 

Total points = 6253; PC= 67.5%; EBD= 2.8% 


Total points = 6253; PC= 77.7% ; EBD= 1.9% 

POD(NE)= 81.7%; POD{SE) 

=65.5% 



POD(N)= 92.8%; POD{E)= 40.1% 



POD(SW)= 27.9%; POP>(AW)=41.8% 


POP>(S)=69.5%; POD{W)= 27.7% 




In Table 1171 the samples are divided in the following 
4 quadrants: (1) wind blowing from the direction between 
0° and 90° (NE winds); (2) wind blowing from the direc¬ 
tion between 90° and 180° (SE winds); (3) wind blowing 
from the direction between 180° and 270° (SW winds); and 
(4) wind blowing from the direction between 270° and 360° 
(NW winds). These quadrants are reported in black in Fig. [3] 


PC is equal to 67.7%. The probability of detection of the 
wind between 0° and 90° (NE winds) is around 82%; be¬ 
tween 90° and 180° (SE winds), it is 65.5%; between 180° 
and 270° (SW winds) it is 27.9%; between 270° and 360° 
(NW winds), it is around 42%. The model performance is 
very good in the quadrant in which the wind speed flows 
more frequently (between 0° and 90° with 82%). The high 
frequence of occurence of the NE winds at Cerro Paranal is 
confirmed by the windrose of Figure [4] (left), between 1997 
and 1999. The worst POD is found for the third quadrant 
(27.9%), but this corresponds to only less than 4% of the 
observed wind total sample at Cerro Paranal. This low fre¬ 
quence of occurence is also confirmed by the windrose of 
Figure 0 (left). In the other quadrants the model perfor¬ 
mance is, in most cases, well above the random case (25%). 
We note that the low frequency in the quadrant [180°, 270°] 
is confirmed by the climatology, and is not limited to our 
sample of 129 nights. If we divide now the sample in 4 dif¬ 
ferent intervals (N, E, S and W, i.e. rotated of 45° with 
respect to the previous one in order to better constrain the 
angular sectors), we observe that the results remains very 
similar. In Table [TH1 this time the samples are divided in the 
following 4 quadrants: (1) wind blowing from the direction 
between -45° and 45° (N winds); (2) wind blowing from the 
direction between 45° and 135° (E winds); (3) wind blowing 
from the direction between 135° and 225° (S winds); and 
(4) wind blowing from the direction between 225° and 315° 
(W winds). These quadrants are reported in red in Fig. 0 
We notice that PC remain high, as expected, with a value of 
77.7%. POD are very good (as high as 92.8% for POD(N)), 
except for POD(W) (less than 28%). In all cases, EBD is 
very low, always inferior to 2.8%, which confirms the good¬ 
ness of the prediction by the model. 

POD(SW) (Tables fTTl) and POD(W ) (Tables IT5T) are the 
only PODs giving values close than the random case, we 
can deduce that the wind blowing between 225° and 270° 
(hatched section of Fig. 0) is reconstructed with the most 
apparent difficulty. However, as previously mentioned, the 
observed frequency of these events is very low. 
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Figure 4. Left panel: observed wind direction at 10 m above the ground at Cerro Paranal. The sectors are in steps of 11.25 degrees, 
the length of the sector represents the fraction of time when the wind comes from this direction. The figure is extracted from the ESO 
astroclimatic webpage at https : / /www.eso.org/gen — fac/pubs / astclim/paranal/wind. Right panel: the probability of detection (POD) 
of the wind direction by the model is summarized by sectors. The circle has been constructed by combining the results from Tables 1171 
and 1181 


The results for the probability of detection of the wind di¬ 
rection by Meso-NH are summarized in Figure [4] (right). 


9 CONCLUSIONS 

In this study, we have quantified the quality of the predic¬ 
tions of the absolute temperature, the wind speed and the 
wind direction in the surface layer 0-30 m at the ESO site 
of Cerro Paranal, by the atmospherical non-hydrostatical 
mesoscale model Meso-NH. First, we identified the bias 
and RMSE bet ween mod el and obser vations. With respect 
to the study of lLascaux et al.l J2013I) , the sample is more 
homogeneous (nights distributed on different periods of the 
year), and richer (129 nights instead of 20). This permitted 
us to obtain very reliable and robust values for the bias, 
the RMSE and the bias-corrected RMSE a. 

To quantify the quality of the model predictions, simula¬ 
tions of 129 nights (distributed between 2007 and 2011) 
have been performed. I 11 addition to this, we also have 
constructed contingency tables and analyzed different 
parameters deduced from these tables, more precisely the 
PC (percent of correct detection), the POD (probability 
of detection of a single event), and the EBD (percent 
of extremely bad detections). These contingency tables 
(and the associated statistical parameters) allowed us to 
quantify the performance of the Meso-NH model above 
Cerro Paranal. 

For the absolute temperature, the percent of correct 
detection (PC) computed from 3x3 contingency tables is 


excellent. It is around 75% range depending on the height 
above the ground (2 m and 30 m). POD are mostly larger 
than 66%. It is and as good as 93.7% at 30 m for the 
detection of temperature inferior to 11.5°C. The worst 
performances are observed in detecting the temperature in 
the narrow range [11.5°,13.5°] at 30 m (POD = 55.2%). 
At every levels, the EBD (extremely bad detections) are 
almost all equal to 0%. In all cases one can note that 
the PC and the POD are well greater than 33% (typical 
value for a random distribution for a 3x3 contingency table). 

For the wind speed, the percent of correct detection 
(PC) computed from 3x3 contingency tables is good 
(around 60%, when the best configurations, AX=100 m, is 
used). The probability to detect the wind velocity in the 
three sub-samples limited by the tertiles of the cumulative 
distribution is good too. The strongest wind are especially 
well detected by the model: POD3 is equal to 74% (at 
10 m) and is equal to 79% (at 30 m). 

In all cases the PC and the POD are greater than 33%, 
typical value for a random distribution. We also proved 
that the configuration with the horizontal resolution 
of 100 m in the innermost domain provides a much bet¬ 
ter model performances in reconstructing strong wind speed. 

For the wind direction, differently from the temperature and 
the wind speed, we have analyzed 4x4 contingency tables. 
If we consider the separation in the four categories (NE, SE, 
SW, NW), the PC is 67.5%. If we consider the separation 
in the four categories (N, E, S, W), the PC is 77.7%. 
Both values are much better than for a random forecast 
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(PC= 25%) of a 4x4 contingency table. The POD of the 
wind is very good in the quadrants from which the wind 
flows more frequently (N and NE) with PODs of the order 
of 92.8% (N) and 81.7% (NE). The POD is not satisfactory 
(~ 27%) only in the narrow angular sector [225°, 270°]. 
However the occurrence of wind blowing from SW is the 
less probable (only 4% of the time for the observed sample 
of 129 nights). The impact of the model performances on 
this case can be considered therefore negligible. 

Concerning the relative humidity, we have shown (in 
the appendix) that the model is able to reconstruct the 
high RH (as high as 80% or more), condition for which the 
dome is closed. 

We therefore conclude that the model performances in 
reconstructing the absolute temperature, the wind speed 
and direction in the surface layer [0,30m] above the ground 
is very satisfactory and results of this study guarantee a 
concrete practical advantage from the implementation of 
an automatic system for the forecasts of these parameters. 
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APPENDIX A: CUMULATIVE 
DISTRIBUTIONS 

We report in this section the cumulative distributions of 
bias, RMSE and a for all parameters and at all levels at 
Cerro Paranal, for the sample pf 129 nights. 



16 F. Lascaux et al 


PARANAL - 129 nights 
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Figure Al. Cumulative distribution of the single nights bias, RMSE and bias-corrected RMSE for the absolute temperature (top row), 
the wind speed (middle row, AX = 500 m and AX = 100 m configurations) and the wind direction (bottom row), at 2 m (temperature) 
and at 10 m (wind speed and direction). 
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Figure A2. Cumulative distribution of the single nights bias, RMSE and bias-corrected RMSE for the absolute temperature (top row), 
the wind speed (middle row, AX = 500 m and AX = 100 m configurations) and the wind direction (bottom row), at 30 m. 
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APPENDIX B: RELATIVE HUMIDITY ing these critical high values, even though the frequency of 

occurrence of such an event is very low. Practically, the as¬ 
tronomer needs to know if the humidity during the night 
will be higher than 80%, a threshold above which the dome 
must be closed, and thus observations are not allowed. Our 
results indicate that the model can give a positive forecast of 
the night relative humidity useful for astronomers. Of course 
it would be of interest to identify a richer sample of nights 
with relative humidity as high as 80% across several years 
and study the behavior of the model on a richer sample. 
This might permit us to have a more robust conclusions on 
this specific subject. 


We report in this section the forecasted temporal evolutions 
for 3 nights in 2011, and 4 nights in 2013, of the relative hu¬ 
midity. These nights are all the nights of 2011 and 2013, in 
which a relative humidity (RH) larger than 80% has been ob¬ 
served. This is the unique interesting condition to be studied 
for astronomical application in this region for these wave¬ 
lengths. As already remembered, in this region the RH is 
mainly of the order of 20-30 % and measurements for such 
a low values are not necessarily reliable. Figure [X3] displays 
the temporal evolution of the relative humidity for these 
nights. The forecasted relative humidity is, in these cases, 
very close to the observed one. In two nights over three in 
2011 the forecasted RH is moreover larger than 80% as well 
as the observed one. In one night (5 July 2011), the fore¬ 
casted RH is slightly inferior 70%) but still very high. 
The four night of 2013 with a high observed RH are well 
reproduced too by the model (predicted values always be¬ 
tween 60% and 100%) even if a small dry bias is present. 
This means that apparently, the model succeeds in predict- 
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From 00 UT to 09 UT - 05 July 2011 
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Figure A3. Temporal evolutions during the night of the relative humidity (in %) at 2 m (left) and 30 m (right) at Cerro Paranal, for 3 
nights in July 2011. In black: the observations; in red: the model outputs. 
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From 00 UT to 09 UT - 06 February 2013 
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Figure A4. Temporal evolutions during the night of the relative humidity (in %) at 2 m (left) and 30 m (right) at Cerro Paranal, for 4 
nights in 2013. In black: the observations; in red: the model outputs. 
















